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The Caldeira-Leggett model describes a microscopic quantum system, represented by a 
harmonic oscillator, in interaction with a heat bath, represented by a large number of 
harmonic oscillators with a range of frequencies. We consider the case when the system 
oscillator starts out in the ground state and then thermalizes due to interactions with the 
heat bath, which is at temperature 6. We calculate the position autocorrelation function 
{x{t')x{t)) of the system oscillator at two different times and study its behaviour in 
the small and large time limits. Our results show that the system oscillator thermalizes 
as expected. We also confirm by explicit calculation that the position autocorrelation 
function exhibits periodicity for imaginary values of the time difference t' — t = ir at 
late (real) times t. 

1. Introduction 

The Langevin equation provides a good description of many dissipative classical 
systems. One approach to finding the quantum analogue of the Langevin equation 
is to quantize the system interacting with an environment as a whole, and then to 
eliminate the environmental degrees of freedom. The interaction between the system 
and the environment gives rise to dissipation. The Caldeira-Leggett (C-L) modeP 
follows this approach, modeling the environment as a large collection of thermal 
harmonic oscillators interacting linearly with the system, which is also modeled as 
a harmonic oscillator. The interaction between the system and the environment is 
assumed to be linear in both coordinates. The parameters of the quantum model 
are set after comparison with the classical equation in the appropriate limit. 

Because the C-L model is simple enough to be analytically solved, it is possible 
to study the approach of the system oscillator to thermal equilibrium in detail and 
to explore how the properties of thermal equilibrium are attained at late times. For 
example, we will ask how the periodicity of Green functions in imaginary time of 
the thermal system emerges from the real-time evolution of the system interacting 
with the environment. We show that the periodicity only holds on the compactified 
interval < t < 1/9, where is the temperature of the environment. For con- 
venience we choose natural units, in which the Boltzmann constant and Planck's 
constant are set to one: H — kg ~ 1. 
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2. The C-L Model 

Consider a particle of mass M in a potential v(x) and in contact with a reservoir 
consisting of large nmnber of harmonic oscillators each with mass m and frequencies 
Wfe. Let X be the position of the particle and Rk be those of the reservoir. The 
interaction between the particle and the environment is assumed to be linear in the 
coordinates. The Hamiltonian of the model is: 

= 1^ + -(x) + X ^ CR, + E ll + E l^-^l^l (1) 

where v{x) = ^muj^x'^ is the potential. The particle trapped in the potential v{x) 
is called the "system" . The reservoir is assumed to be in thermal equilibrium at a 
temperature 9 and acts as a heat bath. 

Using the notation p{x, y, t) for the (x, y)— element of the reduced density matrix 
of the system, the reduced density operator at a time t is given by 

Pix,y,t) = J dx'dy'j{x,y,t,x' ,y' ,0)p{x' ,y' ,0). (2) 

The propagator J of the density operator can be obtained using the Feynman- 
Vernon influence functional method (see eq. (3.38) of Ref. [T]): 

J{x,y,t;x' ,y' ,0) = J DxDyexp{iSo[x,y] - Si[x,y]) (3) 

with 

So[x, y] = Sr[x] - Snly] - M7 / (xx) - yy + xy - yx)dT (4) 



Sb = 
Si[x,y] = 



Jo 

-Mx^ - v{x) + -M{Aojfx^ 

I dijj' J coth — 

V2t' 



dT (5) 



TT 

X / / [x(t) — ?/(t)] cos[w'(r — s)][x(s) — ?/(s)](iTc?s, (6) 
Jo Jo 

where 7 is the relaxation constant of the system, and £7 is the upper frequency 
cutoff for the reservoir. The natural frequency lo of the oscillator is modified by the 
interaction. The renormalized frequency of the system oscillator is 

ujR =u? - (Atj)^ (7) 

with 

(A-)^ = (8) 

TT 

In addition, it is useful to introduce the reduced frequency of oscillation of the 
system in the presence of damping: 

w^ujb" - 72. (9) 
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In general, the relaxation constant 7 depends on the couplings Ck in a compli- 
cated manner. If the environment consists of N oscillators and every C]^ — C, using 
equations (3.34) and ( 3.36) of^ , we get 



7 ^ 



(10) 



3. Approach to thermal equilibrium 
3.1. System behaviour at the large times 

For an undamped harmonic oscillator of frequency w, the propagator of the density 
operator is given by (see eq. (6.26) of Ref. [1]): 

J = F^t) exp [i {K2{t)Xfif + Ki{t)X,^, - L{t)X,if - N{t)XfQ] (11) 



X exp [- (A(i)O' + S(i)Oe« + C{m^))] 



where 



F\t) 



AC ^ 
lim l^(27r)^TT 



(12) 



(13) 



A, B,C, K, L, N are functions of time and K1/2 = K zt Mj/2. Here, J has been 
written in terms of the sum and difference coordinates X = x + y and ^ ~ x — y 
respectively, instead of x and y. 

The initial density matrix for a simple harmonic oscillator in the ground state 

is 



p{xi,y^,0) = 



1 



I-KXq 



■ exp 



(14) 



where xq — (Muj) In terms of the sum and difference coordinates X and ^, 



p(X„C„0) = 



/TTXq 



■ exp 



Axl 



(15) 



The density matrix at time t is given by 



^^eM^K.X.^f-Aer) 



dXid^i exp 



Muj 



X,^,{iKi) 



X exp [X,{~tL^f) + U-iNXf - B^j)] 
This Gaussian integral can be solved by applying the general formula: 

J A 1 J 



exp 



X • A • X 



J X 



cTx 



(27r)" 
det A 



exp 



(16) 



(17) 
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where x and J are vectors and A is a symmetric (n x n) matrix. For n = 2, we 
obtain the density matrix in the final state as 

p{Xf,£,f,t) = FVl67rAfw/0exp [-{a^/ + bX + icX f£,f)] (18) 

where 

z ^ 4:Ki^ + Muj{AC + Muj) 

ACL^ - B^Muj + L^Muj + ABKiL 
z 



b 



c 



AKiLN -2BMNUJ 



(19) 
(20) 

(21) 

(22) 



Thus the density matrix element p{x,y,t) is a Gaussian. 

Plots of the variation of a, b, and c with time are shown in Figs. [T] [2] for the 
parameters 



10 MeV, 7 = 0.2 MeV, n = 100 MeV, M = 130 MeV. 



(23) 



We observe that a, b and c stabilize after some time as what one would expect 
for a system that thermalizes. We define the relaxation time to be ti^i = (27)"^^^, 
which for our choice of parameters is t^ci = 2.5 (MeV)^^. The system is very close 
to thermal equilibrium after about 8 relaxation times, i.e. after about 20 (MeV)~^. 





Fig. 1. The coefficients a{t) and b{t). Botii coefficients attain a constant value at late times. 



Using the normalization condition for the reduced density matrix p{t) , we can 
express F in terms of N as follows. Since the density matrix at t = is normalized 
to one, unitarity dictates that it continues to remain normalized at any time t: 

Tr[p(t)] = / dxp{x,x,t) = 1 (24) 



From ( 18 1 we get: 

Pf{x,x,t) = pj:{X,0,t) = F'^{t)^/l6'KMuj/zexp ^- ^'^^ ^ 



(25) 



December 17, 2012 2:7 WSPC/INSTRUCTION FILE main 



50 



o 



-50 




10 20 30 40 50 60 
t 



70 



Fig. 2. Oscillatory behaviour of the coefBcient c(t), which approaches zero at late times. 



upon substitution into (24) we obtain 



dx ^/l6^^MuJ/z exp 
resulting in the normalization condition 



2^ 



(26) 



(27) 



3.2. Harmonic oscillator at a temperature 9 

Consider now the system oscillator in thermal equilibrium with temperature 9. Its 
density matrix is given by^ 

Mw 



27rsinh(w/6') 
Mw 



exp 



((x^ + y^) cosh(w/6') - 2xy) 



(28) 



2sinh(w/6l) 

In terms of the sum and difference coordinates X and ^ the density matrix can be 
written as 

Mw 



where 



2tt sinh(w/6') 
Mw 



ath = 
hh = 



4sinh(w/6') 

Mw 
4sinh(w/6') 



exp (-flth'C - bthX 



{cosh{w/0) + 1), 
(cosh(w/6i) - 1) 



(29) 



(30) 
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3.3. Comparison of the C-L model with a thermal system 



We now can compare the values of a and b for the thcrmalized system with those 
obtained for the model at late times { t = Sirci = 20 (MeV)^^) . Comparison of 
a(t) and b{t) given by eqs. (20 21) with ath and 6th given by eq. (30) confirms that 



the system oscillator in the C-L model thermalizes, indeed. Table 1 provides a com- 
parison of the dynamically obtained coefficients with the thermal ones for various 
values of 7. Table 2 gives a similar comparison for various values of temperature 9. 



Table 1. Comparison of the C-L model at t = Strci 
with a thermalized oscillator for various values of the 
damping rate 7. 



7 


a 




b 


^'th 


(MeV) 


(MeVf 


(MeV)^ 




(MeV)^ 


0.10 


700.149 


703.285 


151.659 


150.188 


0.15 


708.309 


703.285 


150.452 


150.188 


0.20 


711.141 


703.285 


150.329 


150.188 


0.25 


713.235 


703.285 


150.353 


150.188 


0.30 


715.222 


703.285 


150.403 


150.188 


0.35 


717.187 


703.285 


150.464 


150.188 



Table 2. Comparison of the C-L model at t = Strei 
with a thermalized oscillator for various values of the 
temperature d. 



e 


a 


ath 


b 


feth 


(MeV) 


{MeVf 


(MeV)^ 




(MeV)^ 


10 


711.14 


703.29 


150.33 


150.188 


15 


1015.76 


1010.85 


104.57 


104.492 


20 


1329.25 


1326.97 


79.655 


79.599 


25 


1646.94 


1646.61 


64.192 


64.147 


30 


1966.68 


1968.02 


53.708 


53.670 



4. Position correlation function 
4.1. Two-time correlation function 

We now calculate the two point correlation function of the position operator for the 
system oscillator at different times {x{t + i)x{t)). The time evolution of the density 
operator in the Schrodinger picture can be expressed in two different ways: 

Psit) = Trenv [C/(i)Psys(0)Pcnv(0)C/^(i)] , (31) 
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where U{t) = exp [—iHt] is the time evolution operator and 

Psys{t) = J(t)psys(0). 



(32) 



The density matrices of the system and the environment arc assumed to factorize 
at i = 0. 

We can therefore write the {x, y) element of the density matrix in the following 



two equivalent forms. First, simplifying (31 1, we obtain 



{x\Psys{t)\y) = Trc 



{x \U{t)\ x') (x Vsys (0) Pcnv (0) {v' \ it)\y) dx' dy' 



dx'dy'{x'\psyM\y')T^^cnv [{x\Uit)\x')p,^,{0){y'\UHt)\y)] (33) 



Alternatively, from eq. ( 32 ) we get 



x|peys(t)|y> = / Jix,y,x' ,y' ,t,0){x'\psysmy')dx'dy' 



(34) 



Comparing these two forms, we obtain the following expression for the propagator 
J of the density matrix in terms of the time evolution operator: 



J(x,2/,x',2/',t,0) -Trenv [{x\U{t)\x')p,nAO){y'\U\t)\y)] 



(35) 



We can use this relation to express the time evolution of the Heisenberg picture 
operator x{t) in terms of the propagator J. Since Trenv [penv(O)] = 1, 



{x\x{t)\y) = (a;|x(i)|y)Tro„v [penv(O)] 

= Tronv [{x\x(t)pcnviO)\y)] 



(36) 



Substituting x{t) = W{t)x{0)U{t) in Eqn 36 



we obtain : 



{x\x{t)\y) = Trenv 

— Tr 

— Tr 

— Tr 



dx'dy'{x\U\t)\x'){x'\x{0)\y'){y'\U{t)\y)p,^^{Q) 
dx'dy'{x\UHt)\x')x'S,,^y,{y'\Uit)\y)p,n,{0) 

dx'x'{x\UHt)\x'){x'\U{t)\y)p,nAO) 
dx'x'{x'\U{t)\y)p,^^{0){x\U\t)\x') 



In view of ( 35 1 this implies 



{x\x{t)\y) = / dx'x' J{x' ,x' ,y,x,t,0) 



(37) 



(38) 
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We are now ready to calculate the two point correlation function of the position 
of the system: 

{x{t + i)x{t)) = Tr [p{0)x{t + i)x{t)] 

= Tr,ys [Trenv {piO)U\t + i)x{0)Uit + i)U^ {t)x{0)U{t)}] 
= Tr,ys [Trenv {Uit)p{0)UHt)U^ {i)xiO)Uii)xiO)}] 
= Tr,ys [Tr.nv {p{t}x{i)x{0)}] 



dxdy{x\p{t)\y) {y\x{i)\x)x 



dxdy{y\x{t)\x)x 



Tronv [{x\p{t)\y)] , 



(39) 



where p denotes the density operator for the total system including both, trapped 



particle and heat bath. Using Eqs. (31 38 1, we get 



{x{t + t)x{t)) = / dxdydx' X x' J{x' , x' , x,y,t,0) psys{x,y,t). 



(40) 



Inserting the explicit expressions for J and p, eqs. (12) and (18), we thus obtain 



{x{t + t)x{t)) ^ / dXd^dX 



-p{X,^,t)J{X',0,X,C,t,0) 



= I dXd^dX' ^'^^^^^ exp [- {b{t)X^ + a{t)^^ + tc{t)X^)] 



xM'{t)F'^{i) exp [iKi{i)X^ - iN{i)X'^ - C{i)i^ 



(41) 



where M'{t) = F'^{t)^lQTTMuj / z{t) and a, b and c are given by eqs. (^20 2ll 22L 
The Gaussian integral that can be solved: 



{x{t + t)x{t)) = - 



F^(i)F^{t)y/M^TT^ 



Ki{i)-c{t) i 
2b{t)N{t) N{i) 



(42) 



Using the normalization condition for the density operator we finally obtain the 
autocorrelation function for the position operator: 



{x{t + i)x{t)) = 



2c{t) - Mwcot{wi) - Mj 



where w is the natural frequency of the oscillator and w is given by 

2 2 /^4r2 

w = (jj ^ 7 I ^ 7 



(43) 



(44) 



and and z, b and c are given by eqs. (|l9) [2l] [22]) . 
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4.2. Position uncertainty 



We first consider the equal-time position uncertainty. Substituting i = in eq. (43), 
we find 

{x^{t)) = '-^ (45) 

It is instructive to check the correctness of this resuh at i = and in the hmit 

f — )■ 00. 

For t = 0, {x'^(t)) given by eq. (45 1 has to reduce to the position uncertainty for 
the chosen initial state of the system oscillator. The position uncertainty of a free 
harmonic oscillator in its ground state is given by 



On the other hand, inserting t = in (45), we find: 

1 4:Ki{tf + APuj^ + AMuC{t) 
~ 2Muj 

sm^uit) 

Using C(0) = Q,Ki{t) ^{Muj cot ut + M7) for t ^ above, we get 
1 (Afwcoswt + M7sin^a;t) + APuj^{smujty 



2Muj APlo^ 



1 APuj^{cos^ujt) 



2Muj M^Lu^ 



1 



(47) 



which agrees with ( 46 ) . 

In the limit t — > 00, the position uncertainty of the system oscillator should re- 
duce to that for the harmonic oscillator in thermal equilibrium with its environment. 
For a simple harmonic oscillator in thermal equilibrium, the position uncertainty is 
given by 

2 ^ Tr(a;2pth) ^ smh{w/0) ^ coth [w/{2e)] 

Tr(pth) 2Afu;(cosh(u;/6') - 1) 2Muj ^ ' 



On the other hand, taking the long time limit of (45), we obtain: 

, 2. ^^ 4Xi(i)2 -f 4A/wC-|- Af2a;2 2 

Inn (.^(i)) = 8M..A./2^2 (-2704(sm«;0^ (49) 

Using the explicit forms of K\(t) and C(t) the leading term is: 
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In the limit 7 — >■ the last factor becomes a delta function. Further taking ft 00 
we obtain 

'n I 

'kS{v — w) 



{-m^j^j^ Vdv coll. - 

_ coth [u;/ (26*)] 
~ 2Mw ' 



{v + wY + 7^ 



(51) 



which agrees with ( 48 ) 



4.3. Unequal-time correlator 

We now turn to the unequal time correlator {x{t')x{t)) of the position of the system 
oscillator. Interpreting the quantum oscillator as a (0+l)-dimensional quantum field 
theory, this quantity is the analog of the two-point function of the field operator. 



4.3.1. Long-time limit 

We find the two point correlation function {x{t)x{0))t\\ for a harmonic oscillator in 
thermal equilibrium with a reservoir at a temperature 0: 

(a;(i>(0))th = Tr [p,,-,U\t)x{G)U{i)x{Q)\ 

dxdydz X z pthix,y){y\U^ {i)\z){z\U{i)\x) (52) 

The density matrix of a thermal harmonic oscillator in position space isP 
Pth{x,y) = {x\exp{-/3H)\y) 

Mw 



Mw w 

tanh — exp 

TT 20 



(x^ + y'^) cosh ■ 



2sinh(w/6l) 

The time evolution operator has the position matrix elements: 



2xy 



{z\U{i)\x) 



Mw 



— exp 



iMw 
2 sin wt 



= ((z^ + x'^) coswi- 2xz) 



.(53) 



(54) 



2ni sin wt 

Using these in ( |52| ) and solving the Gaussian integral, we get 

{x{t + i)x{t))th ^ {x{i)x{0))th ^ ^^^j^ cos (wf) coth (^^^ -i sin (55) 

where, as before, w is the effective frequency of oscillation of the damped system 
oscillator. 



We now compare the long-time limit ( 43 ) of the C-L system oscillator with the 

(56) 



thermal oscillator. For the parameter choices 

6* = 10 MeV, 7 = 0.2 McV, fl = 100 MeV, M = 130 MeV, 
we find 



{x{t + i)x{t))th = 0.001736 (MeV)" 
(a;(t + i)x(t))sys = 0.001734 (MeV)" 



(57) 
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for t = -i/b (MeV)"^ and t = 20 (MeV)"\ 

Thus we see that the autocorrelation function of the position operator at differ- 
ent times shows the expected behavior. 

5. Imaginary time correlation function 
5.1. Periodicity of the thermal correlator 

It is well knowrP that the two-point function of any Heisenberg operator O, 
(O(f))O(0)) for a thermal system at a temperature is periodic in imaginary time 
t — —IT with the period /3 — 1/9. Consider the thermal average of the unequal-time 
correlation function of any operator O: 

(O(i))O(0))th = Tr [pthO(i))O(0)] 

= Tr [exp i-PH) exp {iHi)0{Q) exp (-iM)O(O)] (58) 

For t = —i/3 one thus obtains, making use of the invariance under cyclical permu- 
tation of the operators under the trace: 

(O(-i/3))O(0))th = Tr [exp exp (/?i/)O(0) exp (-/3i/)O(0)] 

= Tr [0(0)) exp {-(3H)O{0)] 

= Tr [exp {-f3H)Omom = (O(0)O(0))th. (59) 

We note in passing that this property only implies periodicity of the correlation 
function on the toroidally compactificd interval t G [0, /3], not on the entire imagi- 
nary time axis. For example: 

(O(-2i;9))O(0))th = Tr [exp (-/3-ff) exp (2/3iJ)O(0) exp (-2/3iJ)O(0)] 
= Tr [exp (/3ff)O(0)) exp {-2/3H)O{0)] 
= Tr [exp {~2l3H)O{0)) exp {f3H)O{0)] 
= Tr [exp (-/3iJ) exp {-PH)O{0)) exp (/3i?)O(0)] 
= (O(^/?)O(0))th 

^ (O(-z/3))O(0))th. (60) 

In fact, one easily checks that the "periodicity" property generalizes to the state- 
ment that the thermal imaginary time correlation function is symmetric under re- 
flection at the point t = — i/3/2, rather than to periodicity along the full imaginary 
time axis: 

(O(-z/3/2 + i))O(0))th 

= Tr [exp {-PH) exp (/3H/2 + iHi)O{0)) exp (-/3i7/2 - iHi)O{0)] 

= Tr [exp {-PH/2 + iHi)0{Q)) exp {-PH/2 - iHi)O{0)] 

= Tr [exp i-(3H/2 - iHi)0{Q)) exp (-/3i?/2 + iHt)O{0)\ 

= Tr [exp i-PH) exp (/3i//2 - iHi)0{Q)) exp (-/3i?/2 + iHt)0{Q)\ 

= (0(-z/3/2 - t)O(0))th. (61) 
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5.2. Periodicity of the correlator in the long-time limit 

If the system oscillator thermalizcs in the limit i — > oo, the periodicity in imaginary 
time t must also be valid for the two-time correlation function {x(t + t)x{t)) for the 
system for large values of t: 



lim {x{t - i/3)x{t)) = lim {x{t)x(t)). 



(62) 



We can easily check the approximate periodicity numerically for a large real 
value of t: For example, choosing 

6* = 10 MeV, 7 = 0.2 MeV, Q = 100 MeV, M = 130 MeV, t = 200 (McV)"^ 

we find: 

{x{t)x{t)) = 0.008323 (MeV)"^ 
{x [t - i/e) x{t)) = 0.008323 (MeV)~^ 

An analytical proof can be obtained as follows. Substituting t 
and using the fact that c{t) — !■ as < — >■ cx), we find 



(63) 



{x (t-iP)x{t)) 



A^yb{t)z{t) \smwt\ 



Mw cot {—iwf: 



-iP in (|43|) 
M7" 



Substituting b{t)z{t) = MN^{t)uj in ([21]) we get: 

(x it - iP) xit)) = - — ^ ^ 

Further using the expression N{t) 



I sin wt\ 

Mwe'<^ 



-i + 



Ah{t) 

Mw cot {~iw/3) + M7 



and simplifying, we get: 



{x {t-il3)x{t)) 



2Mw 



- sinh {wP) ■ 



Mw cosh {wP) — iM^ sinh (w/S) 



46(t) 



Now, 



4&(i) AN^Mlu 



^2{x\t)). 



.(64) 



.(65) 



(66) 



(67) 



Using the long-time limit (51 ) for the position variance and taking the limit 7 — > 0, 

1 



(66) becomes: 

(a; {t-iP)x{t)) 



2Mw 

coth iwP/2) 



2Mw 
coth{wP/2) 
2Mw 



[coth {wP/2) cosh {wP) — sinh {wP)] 

cosh-{wP) — tanh(w/3/2) sinh(?i;/3)] 



(68) 



In the last step we used the relations cosh(a;) = 1 + 2sinh'^(x/2) and sinh(a::) = 
2 sinh(a;/2) cosh(a;/2). This coincides with the result (51) for (a;^(t)) and thus con- 
firms the periodicity of the two-time correlation function in the long-time limit. 
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5.3. Correlation function in imaginary time 

The Figures |3] - 14] below explore the behavior of the two-time correlation function 

G{t,i) ^ {x{t + i)x{t)) (69) 

in the Caldeira-Leggett model for imaginary values of the time difference t = —it. 
The differently colored lines represent the system oscillator at various times t. The 
black dashed line always represents the thermal system. 



0.0008 
-p 0.0006 




- 0.0004 




0.00 0.02 0.04 0.06 0.08 0.10 

T 



Fig. 3. Periodicity of the two-point function G{t, —ir) on the interval < T < /3 for d = 10 MeV. 
The brown, blue, green, and red lines correspond to times t = 0,5,10,20 (MeV)~^, respectively. 
For times of the order of the thermalization time, the two-point function approaches periodicity. 
The thermal correlation function is shown as a black dashed line. 



The approach of the imaginary-time correlation function of the C-L model to- 
wards the thermal correlator is obvious from these figures. The periodicity of the 
late-time two-point function on the compactified interval < r < /3 is clearly visible 
in Fig |3] The smooth approach to the thermal correlator is shown in Fig. |4] 

Figure [5] which shows the imaginary time correlation function in comparison 
with the thermal correlator over a larger imaginary time domain, illustrates the 
absence of periodicity of the correlation function beyond the interval < r < /3. The 
exponential growth of the correlation function for t > P can easily be understood 
by realizing that the thermal correlator can be written in terms of the eigenstates 
\n) of the harmonic oscillator as 

(x(-zr)a;(0))th = ^ \{n\x\m) p exp[-£;„/3 -f {E^ ~ E^)t\. (70) 

m.n 

Since the position operator connects states to = n± 1, this implies that the thermal 
correlator grows like exp(wr) at large imaginary times. 
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Fig. 5. Comparison of the two-point function G{t, —ir) with the thermal equilibrium correlator 
on the extended imaginary time interval < r < 0.4. Near the thermalization time, the two-point 
function for the system merges with that for a thermal system. The absence of periodicity over 
the larger time interval is clearly seen. 



5.4. Correlation function in real time 

For completeness, we also investigate the dependence of the two-time correlation 
function on real time difference. The figures below show the behavior of the two- 
point correlation function {x{t + i)x{t)) for the Caldeira-Leggett model in real 
relative time i. Figure [g] shows the correlation function G{t, r) as function of the 
real time difference r for several values of t. The correlator is seen to approach the 
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thermal value rapidly. Figure [7] shows the difference between the real-time corre- 
lation function of the C-L model and the thermal correlator on a semilogarithmic 
scale. Finally, Fig. [8] shows the real-time correlation function over a more extended 
time interval. 
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Fig. 6. Correlation function G(t,r) as function of the real time difference r in the interval < 
T < f3 for several values of t. The brown, blue, green, and red lines correspond to times t = 
0,5,10,20 (Mev)"\ respectively. 



0.1 

I 0.001 _ 

"P [ 
+ 

O 10 




0.10 



Fig. 7. Deviation of the correlation function for C-L model from the thermal system in real time 
T in a semilogarithmic plot. 
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Fig. 8. Comparison of the two point function for C-L system witli the thermal system in real time 
. Near the thermalization time, the curve for the system approaches that for a thermal system. 

5.5. Two point correlation function for complex time 

For completeness, we present contour plots showing the behavior of the real and 
imaginary parts of the two point correlation function for time t in the complex 
plane in Fig. [9] 




0.0 0.1 0.2 0.3 0.4 0.0 0.1 0.2 0.3 0.4 



Re[T] Re[T] 

Fig. 9. Contour plot of the real part (left) and imaginary part (right) of the correlation function 
G(t, t) as a function of time t in the complex plane. Here, t = 20 (MeV)~^, i.e. the system is 
thermalized. Lighter shading indicates higher values. 
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6. Summary 

We have studied the time evolution of the density matrix of a harmonic oscillator 
in interaction with a thermal environment, starting in the ground state. Our study 
makes use of the explicit solution obtained by Caldeira and Leggett when the en- 
vironment is modeled as a large collection of harmonic oscillators. We confirmed 
that the parameters of the solution approach constant values at late times. We also 
showed that these asymptotic values agree with those for the density operator of 
the harmonic oscillator in thermal equilibrium. Our comparison with the thermal 
system shows that the system oscillator of the Caldeira-Legett model thermalizes 
as expected. 

We then computed the two-time correlation function of the position operator 
for the Caldeira- Leggett model G{t,i) = {x{t + i)x{t), where i denotes the time 
difference and compared it with the position correlation function for a harmonic 
oscillator in thermal equilibrium. We studied the following cases in detail: 

(1) t ~ 0, corresponding to the position uncertainty of the harmonic oscillator rep- 
resenting the "system" . Here we showed that the uncertainty changes smoothly 
from that of the harmonic oscillator in its ground state at i = to that of a 
thermal oscillator in the long-time limit. 

(2) Imaginary t. We found that the correlation function G{t,t) approaches the 
thermal correlation function {x{t)x{0))th m the long-time limit i — > oo and 
exhibits periodicity in imaginary time t, i. e. limt_>oo [G(t, —i/3)/G(t,0)] — 1. 

(3) Real t. Again, we found that the correlation function G{t, t) approaches the 
thermal correlation function {x{t)x{0))th m the long-time limit t — ^ oo, but the 
function does not exhibit periodicity properties along the real t-axis. 

Our results suggest that the thermalization of a subsystem of a larger quantum 
system evolving in real time can be explored by calculating two-time correlation 
functions of the form 

C{t,T)^{0{t-tr)Oit)) (71) 

of operators defined on the subsystem for real values of t and r and then studying the 
behavior of C as a function of r for large values of the real time t. If the subsystem 
thermalizes, C will become periodic in the variable r, i. e. C{t,/3) C{t,0) for 
t ^ oo. The period /3 is a measure of the temperature 9 = 1//3 to which the 
subsystem thermalizes. 

Our results also raise some interesting questions. For example, if one considers 
a large system in isolation that "thermalizes" under its internal dynamics, will 



correlation functions of the form ( 71 1 exhibit periodicity in t for large values of tl 
Will this only be true for certain operators O, which are sensitive only to a small 
subset of all degree of freedom of the large system, e. g. single-particle observables in 
a many-body quantum system? For a small quantum system, can the average over 
many degrees of freedom be replaced with an ergodic average over the observation 
time tl 
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